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Preface 


This  contract  report  summarizes  a  literature  review  of  suspended-sediment 
transport  theories  under  oscillatory  flow  conditions  performed  by  Dr.  Thomas 
G.  Drake.  The  work  was  performed  under  Evans-Hamilton,  Inc.,  Contract  No. 
DACW39-92-D-0012  with  the  U.S.  Army  Engineer  Waterways  Experiment 
Station’s  (WES)  Coastal  Engineering  Research  Center  (CERC). 

This  report  supplies  a  list  of  suspended-sediment  transport  theories  to  be 
tested  with  field  data  collected  by  personnel  in  the  Prototype  Measurement  and 
Analysis  Branch,  CERC.  Dr.  Thomas  White  was  Principal  Investigator  and 
provided  technical  oversight  for  Work  Unit  32685  "Field  Tests  of  Sediment 
Transport  Theories"  in  the  Coastal  Research  Program  at  CERC.  The  work  unit 
continues  under  the  Inlets  Research  Program. 

WES  supervisors  during  this  report’s  preparation  were  Branch  Chief 
William  L.  Preslan,  Division  Chief  Thomas  W.  Richardson,  Program 
Manager  Carolyn  M.  Holmes,  Assistant  Director  of  CERC  Charles  C. 

Calhoun,  Jr.,  Director  of  CERC  Dr.  James  R.  Houston,  Director  of  WES 
Dr.  Robert  W.  Whalin,  and  Commander  of  WES  COL  Bruce  K.  Howard, 
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1  Introduction 


Background 

Numerous  theories  are  available  to  predict  the  transport  of  cohesionless 
sediments  by  suspension  in  the  surf  zone.  The  purpose  of  this  report  is  to 
survey  the  available  literature  for  such  theories;  and,  subject  to  certain  criteria 
discussed  below,  to  differentiate  between  unique  theories  and  essentially 
redundant  ones.  For  each  of  the  unique  theories,  a  common  set  of  variables  is 
used  to  formulate  an  explicit  mathematical  expression  for  the  suspended 
sediment  transport  rate.  Redundant  theories  are  those  which  differ  only  by 
numerical  coefficients  from  any  of  the  unique  theories,  having  otherwise  the 
same  functional  form.  Accompanying  each  unique  formula  is  an  outline  of  its 
derivation,  and  a  list  of  the  assumptions  underlying  the  derivation.  The 
selected  formulae  will  form  a  basis  for  testing  theoretical  predictions  with  field 
measurements  of  suspended-sediment  transport  and  requisite  fluid  and  sediment 
quantities  at  a  point  within  the  surf  zone  or  immediately  adjacent  to  it. 


Organization 

The  organization  of  the  report  is  as  follows:  first,  the  criteria  used  to  select 
theories  from  the  literature  are  presented;  second,  a  section  on  theory  discusses 
general  considerations,  and  presents  each  of  the  unique  formulas,  their 
derivation,  and  assumptions.  A  table  of  selected  theories  summarizes  formulae 
for  testing  field  measurements.  An  extensive  bibliography  is  presented  in 
addition  to  the  section  describing  references  cited  in  the  body  of  the  report. 
Finally,  Appendix  A  describes  the  notation  used  in  this  report. 


Selection  Criteria 

This  report  emphasizes  theories  for  suspended  sediment  transport  in  or 
immediately  adjacent  to  the  surf  zone  in  combined  wave-current  nearshore 
flows  under  conditions  such  that  quasi-unidirectional  currents  do  not  dominate 
the  combined  flow. 
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The  criteria  used  to  select  theories  for  further  consideration  are  as 
follows: 


a.  Theories  must  be  valid  within  the  surf  zone. 

b.  Theories  must  specifically  predict  suspended-sediment  transport  at  a 
point;  theories  yielding  only  predictions  of  transport  averaged  over 
any  spatial  coordinate  other  than  the  vertical  are  excluded. 

Acceptable  theories  may  thus  require  as  input  depth-averaged  fluid 
and  sediment  quantities,  but  theories  predicting  only  transport 
averaged  over  the  width  of  the  surf  zone,  for  example,  are  excluded. 

c.  Specifically  excluded  are  formulations  using  only  bed  elevation 
measurements  as  input  and  the  kinematic  constraints  of  the  sediment 
conservation  equation  to  calculate  spatial  gradients  in  the  sediment 
transport  rate.  Such  formulations  do  not  provide  predictions  of  the 
sediment  transport  rate  based  on  measured  fluid-motion  and  grain-size 
parameters. 

The  comprehensive  literature  search  focused  primarily  on  theories  presented 
in  refereed  journals  and  government  reports,  because  they  are  readily  available 
in  most  libraries  and  have  been  peer  reviewed.  However,  a  considerable 
sediment-transport  literature  exists  outside  this  group,  primarily  in  the  form  of 
conference  proceedings,  and  these  reports  are  included  in  the  survey  to  the 
extent  that  such  reports  are  available.  Fortunately,  there  is  considerable 
overlap  between  the  refereed  and  so-called  "gray"  literature,  so  that  essentially 
all  useful  transport  formulae  published  in  the  English  language  (or  translated  to 
English  from  the  original  literature)  are  included  here. 
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2  Theory 


Theories  for  transport  of  suspended  sediment  fall  into  three  broad 
categories:  those  based  on  explicit  energy  arguments;  those  based  on  the 
diffusion  theory  or  modifications  to  it;  and  quasi-empirical  predictors  of 
sediment  transport.  With  few  exceptions,  theories  for  suspended  sediment 
transport  in  the  surf  zone  are  directly  based  on  or  are  fundamentally  related  to 
theories  first  developed  to  describe  steady  unidirectional  flows  in  fixed-bed 
channels.  Thus,  the  original  theories  for  unidirectional  flow  are  occasionally 
referenced  to  clarify  the  derivations  of  surf  zone  transport  theories. 

Because  the  so-called  "energetics"  models  require  the  measurement  of  far 
fewer  quantities  to  predict  transport  in  the  surf  zone,  they  are  greatly  favored 
in  nearshore  field  studies.  The  following  sub-sections  outline  the  energy-based 
approach  of  Bagnold  (1963;  1966),  critically  examine  assumptions  made  in  a 
number  of  energetics-based  theories  that  extend  Bagnold’s  approach,  and  list 
the  limitations  of  those  theories;  then  a  similar  approach  is  taken  with  the 
diffusion  theories.  Finally,  two  commonly  used  quasi-empirical  predictors  of 
sediment  transport  are  critically  examined.  Particular  emphasis  is  given  to 
calibration  constants  and  free  parameters  in  the  theories. 


Energetics-Based  Theories  for 
Suspended-Sediment  Transport 


Energetics-based  theories  for  suspended  sediment  transport  following  the 
assumptions  of  Bagnold  (1963;  1966)  have  been  proposed  by  numerous 
investigators.  The  form  of  the  general  theory  and  its  assumptions  are 
presented  below,  followed  by  a  critical  examination  of  the  modifications  to  it 
for  unsteady  oscillatory  flow  proposed  by  Bowen  (1980),  Bailard  (1981;  1982; 
1987)  and  Roelvink  and  Stive  (1989). 

Bagnold  (1966)  derived  relations  for  the  transport  of  sediment  as  bed  load 
and  suspended  load  in  unidirectional  streams  based  on  the  assumption  that  a 
fixed  fraction  of  the  total  stream  power  is  expended  moving  bed  load,  and  that 
a  fixed  fraction  of  the  remaining  available  power  is  expended  transporting 
suspended  sediment.  Bagnold  (1966,  p.  1-4)  defined  bed  load  as  "that  part  of 
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the  load  which  is  supported  wholly  by  a  solid-transmitted  stress,"  and 
suspended  load  as  "that  part  which  is  supported  by  a  fluid-transmitted  stress. 
Taking  the  total  stream  power  0  as 


©  =  p  ghSu  D' 

where 

p  =  fluid  density 
g  -  gravitational  acceleration 
h  =  flow  depth 
S  =  energy  slope 
u  =  mean  flow  velocity 

Bagnold  assumed  that  the  immersed-weight  transport  rate  of  suspended  load  is 


esu(l-eb) 

W 


0 


(2) 


where 

es  =  suspended-load  transport  efficiencies 
eb  =  bed  load 
W  =  sediment  fall  velocity1 


Unfortunately,  there  is  an  error  in  Bagnold’ s  derivation  of  the  suspended 
sediment  efficiency  es  that  appears  to  seriously  compromise  the  theory.  In 
Bagnold’s  Equation  14  (p.  1-14),  he  mistakenly  suggests  that  the  momentum 
flux  due  to  downward-directed  momentum  moving  downward  is  opposite  that 
due  to  upward-directed  momentum  moving  upward,  when  in  fact  both  have  the 
same  sign. 

Despite  its  deficiencies,  Bagnold  and  others  have  extended  the  energetics 
approach  to  predict  sediment  transport  in  unsteady  oscillatory  flows.  For 
oscillatory  flows,  Bagnold  (1963)  proposed  a  similar  relation  for  the  total  load: 


i  =  K  a>  -1 
6 


(3) 


1  For  convenience,  symbols  are  listed  in  the  notation  (Appendix  A). 
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where 


Kb  =  constant 

m8  =  current  velocity  in  the  0  direction 

um  -  maximum  wave  orbital  velocity  measured  at  the  same  distance  above 
the  bed  as  u6 

z'e  =  transport  in  the  0  direction 


Although  not  considered  further  here,  functionally  similar  expressions  for  the 
total  load  have  been  developed  by  many  investigators.  The  following  para¬ 
graphs  outline  and  contrast  the  modifications  to  Bagnold’s  ideas  proposed  by 
Bowen  (1980),  Bailard  (1981;  1982;  1987),  and  Roelvink  and  Stive  (1989). 

Common  to  each  of  the  energetics-based  theories  are  the  following 
assumptions: 

a.  The  immersed  sediment  transport  rate  i  is  proportional  to  to,  the  local 
rate  of  energy  dissipation.  Assuming  that  the  shear  stress  x  is 
quadratic  in  the  velocity  u 


X  =  pcfu 2  (4) 

where  cf  is  the  bed  drag  coefficient,  and  the  rate  of  energy  dissipation  to 
is  assumed  to  be 


co  =  pcy|  u  |  ^ 

b.  Equations  4  and  5  assume  zero  phase  difference  between  stress  and 
velocity. 

c.  The  bathymetry  is  free  of  bed  forms,  and  has  parallel  contours. 

d.  Flow  accelerations  play  no  role  in  the  sediment  transport  processes. 

e.  The  sediment  grain-size  distribution  is  completely  characterized  by  a 
single  constant  fall  velocity  W,  which  is  independent  of  concentration. 

/.  There  is  no  provision  for  initiation  of  sediment  motion. 

g.  Sediment  transport  responds  instantaneously  to  changes  in  flow 
conditions. 

h.  An  implicit  assumption  is  that  the  concentration  of  suspended  sediment 
is  sufficiently  low  that  stratification  effects  can  be  neglected  (e.g., 
Villaret  and  Trowbridge  (1991)). 
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i.  Vertical  fluid  velocities  are  assumed  to  play  no  role  in  the 
depth-integrated  transport  process,  beyond  their  obvious  role  in 
suspending  sediment.  There  is  mounting  evidence  (e.g.  Bowen  (1990), 
summary  in  Nielsen  (1992)  that  vertical  velocities  can  be  important, 
particularly  if  sediment  concentration  decreases  slowly  with  distance 
from  the  bed. 

Generalizing  the  expression  developed  by  Bailard  (1981),  the  time-averaged 
immersed  suspended-sediment  transport  rate  at  a  point  is 


{  I*)  =  pcy-fi [<|  u\3it)  +  Kl+  £2(|m|5>  i] 

where  angle  brackets  indicate  time  averaging,  u  is  the  instantaneous  fluid 
velocity  vector,  f  is  the  unit  vector  in  the  downslope  direction,  and  the 
constants  K}  and  K2  are  discussed  below.  For  purposes  of  discussing  proposed 
variations  of  Equation  6,  the  primitive  equation  (e.g.,  Guza  and  Thornton 
(1985))  for  suspended-sediment  transport  in  the  cross-shore  (x)  direction  in  the 
case  Kx= 0  is 


M  _  escf  p“3N 
W(\-K2u) 

Here  u  is  the  velocity  component  parallel  to  the  local  bed  slope. 


(7) 


Bowen  (1980) 

Bowen  adopted  Bagnold’s  model  as  a  simple  starting  point  for  developing 
models  of  onshore-offshore  transport.  Explicitly  supposing  the  maximum 
orbital  wave  velocity  um  to  dominate  the  total  velocity 

u  =  um  cos  (at)  +  Ul  (8) 

where  a  is  the  radial  wave  frequency,  Bowen  expanded  Equation  8  and 
investigated  several  cases  for  various  forms  of  the  perturbation  Uv 

Both  Bagnold  and  Bowen  derived  £,=0  and  K2=fi/W  (implicitly  assuming 
p«l,  so  that  tanp  =  p,  where  P  =  angle  of  bed  slope).  Bowen  specifies 
temporal  averaging  over  many  wave  periods. 


Bailard  (1981;  1982;  1987) 

Bailard  derived  a  slightly  different  form  of  Equation  6,  keeping  only  the 
first  two  terms  in  the  binomial  expansion  of  Equation  7.  Bailard’s  derivation 
gives  J5T,=0  and  K2  =  e,tanp/Vf.  As  Bailard  (1981,  p.  10,952)  notes,  "this 
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difference  arises  because  Bagnold  assumed  that  the  stream  power  contribution 
from  the  suspended  sediment  load  contributes  directly  to  the  suspended 
sediment  transport  rate  instead  of  through  an  efficiency  factor..." 

The  possibility  of  "autosuspension"  arises  when  the  factor  ( \-K2u )  in  the 
denominator  of  Equation  7  approaches  zero.  The  expansion  leading  to  (6)  is 
singular  in  this  case.  First  noted  by  Knapp  (1938),  autosuspension  refers  to 
turbulent  flows  transporting  suspended  sediment  in  which  the  rate  of  energy 
input  to  the  flow  by  the  downslope  motion  of  the  sediment  equals  or  exceeds 
the  rate  at  which  energy  is  expended  by  the  flow  to  suspend  the  particles.  In 
such  cases,  additional  sediment  can  be  suspended  by  the  flow  with  no  net 
energy  input,  and  thus  the  concentration  is  limited  only  by  sediment 
availability. 

The  models  of  Bailard  and  Bowen  both  admit  the  possibility  of 
autosuspension,  though  the  quantitative  criteria  differ.  Bailard  (1981)  and 
Pantin  (1979)  predict  autosuspension  when 

tanp  >  —  (9) 

esu 

and  Bagnold  (1962)  and  Bowen  (1980)  predict  autosuspension  when 

um|)>-  (10) 

u 

Clearly  such  conditions  can  arise  in  the  surf  zone,  and  the  viability  of  the 
autosuspension  concept  has  been  debated  in  the  literature.  Critical  laboratory 
experiments  by  Southard  and  Mackintosh  (1981)  show  that  the  Bagnold  auto¬ 
suspension  criteria  (Equation  10)  is  not  valid.  Their  experiments,  which  do  not 
test  the  more  conservative  criterion  (Equation  9)  proposed  by  Bailard  and 
Pantin,  strongly  suggest  that  the  energy  balance  in  these  flows  is  not  merely  a 
matter  of  comparing  the  energy  input  by  the  downslope  motion  of  suspended 
particles  and  the  energy  expended  to  suspend  the  particles.  That  comparison 
ignores  the  changes  induced  by  the  presence  of  suspended  particles  on  the 
turbulent  structure  of  the  flow;  and  perhaps  of  greater  importance,  does  not 
include  the  effects  of  high  suspended-sediment  concentration  on  the  boundary 
resistance.  Turbulence  is  damped  by  density  stratification  at  high 
concentrations,  thus  limiting  the  vertical  transport  of  mass  and  momentum 
(McLean  1991).  Thus  boundary  resistance  can  be  significantly  decreased  as 
the  suspended  load  increases,  in  turn  reducing  the  ability  of  the  flow  to  entrain 
more  sediment  into  suspension.  It  is  this  balance  between  entrainment  and 
deposition  at  the  bed,  Southard  and  Mackintosh  argue,  that  governs  the 
equilibrium  concentration  of  sediment  in  the  flow.  Further  discussion  of  the 
issues  can  be  found  in  Parker  (1982),  Pantin  (1982),  Paola  and  Southard 
(1983),  and  Seymour  (1986). 
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Roelvink  and  Stive  (1989) 


The  theory  proposed  by  Roelvink  and  Stive  (RS)  assumes  K2  =  0,  and 
derives  Kx  based  on  a  one-equation  model  for  turbulence  produced  by  breaking 
waves.  The  overall  goal  of  their  theory  is  to  incorporate  models  for  a  variety 
of  cross-shore  flows,  such  as  asymmetric  oscillatory  flow,  wave 
grouping-induced  long-wave  flow,  breaking-induced  turbulent  flow,  and 
momentum  decay-induced  undertow.  The  resulting  model  for  the  fluid  flow, 
derived  for  cross-shore  transport  only,  can  be  simplified  and  used  for 
comparison  with  field  data.  In  their  model,  the  parameter  Kx  is 

tfi=PdUjexp(/l/tfnJ-l]-1}3'2  (ID 

where  pd  is  a  dimensionless  coefficient  (assumed  to  be  =1),  =  root-mean- 

square  wave  height,  and  k ,  is  the  turbulent  kinetic  energy  per  unit  mass 

k  = -L  [(m'2)  +  (v'2)  +  (w'2)]  (12) 

'2 

where  primes  indicate  fluctuating  velocity  components.  Although  the  RS 
model  calculates  k,  by  an  iterative  solution  of  a  system  of  equations,  available 
field  data  may  be  used  to  determine  k,  directly.  Alternatively,  the  average 
dissipation  Ed  per  unit  mass  due  to  wave  breaking  might  be  approximated 
following  Thornton  and  Guza  (1986): 

e,=  —yftTgJL.  /Hi  (13) 

16  y*h5 

where /p  is  the  peak  frequency,  y  is  an  empirically  determined  ratio  of 
root-mean-square  wave  height  H ^  to  the  water  depth  h  in  the  surf  zone 
(Thornton  and  Guza  1983),  and  is  typically  about  0.4  ±  0.1.  The  coefficient  B 
accounts  for  differences  in  wave  breaking;  it  is  less  than  one  for  spilling 
breakers  and  near  unity  for  plunging  breakers  (Thornton  and  Guza  1986). 

Then,  assuming  a  balance  between  the  dissipation  ed  and  the  creation  of 
turbulent  kinetic  energy  gives 


which  supplies  all  the  information  necessary  to  determine  Kx. 


Humiston  (1993) 

Following  Guza  and  Thornton  (1985),  Humiston  modeled  contributions  of 
infragravity  and  incident  waves  and  quasi-steady  currents  to  cross-shore 
transport  of  both  bed  load  and  suspended  load.  The  total  velocity  is 
decomposed  into  quasi-steady  and  oscillatory  components 
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(15) 


ut  =  (u  +  u)i  +  (v+  v)j 

where  overbars  indicate  the  quasi-steady  quantities,  tildes  indicate  oscillatory 
ones,  and  i  and  j  are  unit  vectors  in  the  on-offshore  and  alongshore  directions. 

The  oscillatory  terms  are  further  expanded  into  infragravity  u,  (f<0.05  Hz) 
and  incident  us  (f>0.05)  components  such  that 

u  =  u, +  us  (16) 

and 

V  =  V,  +  VS  (17) 

The  magnitude  of  the  total  velocity  is  then 

|w,|  =  [n2  +  v2  +  u2  +  v2  +  2(u  u  +  v  v) ]7 

This  decomposition  of  the  oscillatory  flow  velocities  may  be  useful  in 
comparing  model  outputs  with  field  measurements. 

Bed  friction  and  suspended  sediment  efficiency  factors 

The  bed  friction  and  suspended  sediment  efficiency  factors  are  two 
important  parameters  in  the  energetics-based  theories.  Although  both  are 
assumed  constant  in  the  theories  presented  above,  it  is  likely  that  the  bed 
friction  factor  is  a  complex  function  of  fluid-motion  parameters,  in  particular, 
the  properties  of  the  bed-load  layer  and  any  small-scale  bed  roughness 
generated  by  the  bed-load  process.  The  choice  of  suspended  sediment 
efficiency  factor  lacks  theoretical  guidance  to  indicate  its  value,  and  must  be 
determined  by  calibration.  The  following  discussion  focuses  on  Bailard’s 
studies,  which  typically  form  the  basis  for  the  choice  of  constants  in  other 
investigations. 

Bailard  (1981)  determined  the  suspended  sediment  efficiency  factor  using 
selected  data  from  laboratory  and  field  experiments.  He  determined  e3  =  0.025 
(as  well  as  the  bed-load  efficiency  factor,  the  bed  friction  factor,  and  a  number 
of  other  parameters)  from  measured  estimates  of  total  longshore  transport  and 
velocity  using  an  iterative  scheme  from  Ostendorf  and  Madsen  (1979). 
Estimated  values  of  the  parameters  are  not  independent,  and  error  analysis  is 
difficult.  Bailard  (1987,  p.  339)  used  es  =  0.035  for  field  data  from  Leadbetter 
and  Torrey  Pines  Beaches  as  part  of  the  nearshore  sediment  transport  study. 
Other  investigators  have  adopted  Bagnold’s  (1966)  estimate  es  =  0.015  for 
unidirectional  stream  flow.  Since  the  transport  rate  is  directly  proportional  to 
es,  errors  of  100  percent  or  more  may  be  expected  in  uncalibrated  application 
of  Equation  6  to  field  measurements  due  to  errors  in  estimating  es  alone. 
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The  measured  values  of  cf  in  the  experiments  Bailard  (1981)  used  to 
determine  es  range  from  0.013  to  0.017  for  the  laboratory  experiments;  for  the 
three  Silver  Strand  field  experiments  (see  Komar  (1969)  for  original  data)  used 
in  the  analysis,  pranged  from  0.0026  to  0.014.  It  should  be  noted  that  the 
values  for  cf  presented  in  Bailard’s  Table  2  (p.  10,951)  are  inconsistent  with 
his  Equations  42  and  43  (p.  10,946)  using  T  =  0.13  and  values  for  other 
parameters  from  his  Table  2.  There  is  no  discussion  of  the  value  adopted  for 
cf  in  the  calculations  in  the  paper,  although  subsequent  papers  (Bailard  (1982), 
p.  1424)  cite  a  value  of  cf=  0.005  based  on  Komar’s  (1969)  Silver  Strand 
data. 

Because  the  bed  friction  factor  cf  is  an  integral  part  of  many  longshore 
current  models,  there  has  been  considerably  more  study  of  cf  than  es.  Shemdin 
et  al.  (1978)  compiled  values  of  c{  for  various  field  studies  of  wave  dissipation, 
and  found  a  range  of  0.005  to  0.3.  Thornton  and  Guza  (1986)  determined  ca¬ 
using  different  longshore  models;  for  a  linear  bed  shear  stress  model  (Equation 
14,  p.  1167)  they  found  the  optimal  cy=  0.009+10.001,  while  for  the  nonlinear 
model  (Equations  22  and  23,  p.  1168)  the  best-fit  cf  =  0.006110.0007, 
significantly  less  than  for  the  linear  model.  Their  linear  model  corresponds  to 
the  quadratic  shear  stress  model  Equation  4  assumed  in  the  energetics-based 
transport  models.  Thus,  uncalibrated  application  of  Equation  6  and  its  variants 
may  introduce  uncertainties  of  100  percent  and  more  into  transport  estimates 
due  to  errors  in  cf . 


Diffusion-Based  Theories  for 
Suspended-Sediment  Transport 


Several  theoretical  approaches  to  suspended  sediment  transport  seek  to 
calculate  transport  based  on  a  detailed  picture  of  near-bed  fluid  motion,  in 
particular,  the  computation  of  the  boundary  shear  stress  and  the  assumption 
that  fluid  turbulence  diffuses  sediment  according  to  a  Fickian  diffusion  law. 
The  assumption  of  Fickian  diffusion  is  not  the  only  one  possible,  but 
alternatives  to  it  (e.g.,  McTigue  (1983))  have  not  been  used  in  the  context  of 
nearshore  transport  and  thus  are  not  considered  further  here. 

Rouse  (1937)  derived  a  diffusion-abased  theory  for  suspended  sediment 
transport  in  steady,  uniform  unidirectional  flow,  which  forms  the  basis  for 
theories  describing  transport  in  unsteady,  combined  flows  in  the  nearshore. 
The  governing  convection-diffusion  equation  for  the  volume  concentration  of 
sediment  N  is 


M  M 

ST  ST 


+ 
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Tz 
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where  us  is  the  sediment  diffusivity,  often  equated  to  the  fluid  eddy  viscosity, 
and  diffusion  and  advection  in  the  x  and  y  directions  have  been  neglected.  For 
steady  uniform  flow  in  a  channel,  the  vertical  distribution  of  suspended 
sediment  is  then 


N=N0 


h-z  zo 


The  Rouse  number  R  is 


(20) 


where  cr  is  an  empirically  determined  constant  of  proportionality  between  the 
sediment  diffusivity  and  the  fluid  eddy  viscosity  (often  assumed  to  be  unity), 
k  is  von  Karman’s  constant,  and  u,  =  (zb  /p)1/2  is  the  friction  velocity.  The 
derivation  of  Equation  19  assumes  a  law-of-the-wall  velocity  profile; 
application  of  it  requires  specification  of  a  boundary  condition  at  z  =  z0. 

Typical  extensions  to  the  Rouse  formulation  for  unsteady  flows  include  a 
time-varying  and/or  vertically  varying  eddy  viscosity,  and  a  bed-shear  stress 
relation  for  the  reference  concentration  N0. 

Kennedy  and  Locher  (1972)  and  Smith  (1977)  present  thorough  discussions 
of  the  diffusion-based  theories  and  their  application  to  both  oscillatory  and 
unidirectional  flows.  Grant  and  Madsen  (1979)  use  an  eddy-diffusivity  closure 
scheme  to  couple  equations  for  conservation  of  fluid  momentum  and  sediment 
mass  to  describe  combined  wave-current  flows  on  the  continental  shelf.  They 
solve  the  model  equations  using  an  iterative  numerical  technique.  Glenn  and 
Grant  (1987)  extended  the  work  of  Grant  and  Madsen  (1979)  to  include 
stratification  effects.  Wikramanayake  and  Madsen  (1992)  used  new  relations 
for  the  reference  concentration  at  the  bed  in  an  otherwise  similar  approach  to 
suspended  sediment  transport  under  unbroken  waves. 

Diffusion-based  theories  have  been  extended  to  treat  several  special  cases  of 
sediment  transport  in  the  surf  zone,  but  there  have  been  no  theories  developed 
to  treat  the  general  case  of  combined  wave-current  flow  that  have  been 
reduced  to  a  single  formula  for  comparison  with  field  measurements;  the 
analyses  typically  result  in  systems  of  coupled  nonlinear  equations  that  must  be 
solved  numerically.  Major  difficulties  in  applying  such  theories  to  field  data 
include  specification  of  initial  and  boundary  conditions  with  sufficient  detail  to 
ensure  unique  solutions  to  the  equations.  A  brief  review  of  diffusion-based 
models  for  suspended  sediment  transport  in  the  surf  zone  is  presented  below; 
however,  none  of  the  theories  meet  the  selection  criteria  and  thus  are  not 
included  in  Table  1. 

In  general,  the  instantaneous  rate  of  suspended  sediment  transport  is  given 
by 
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where  u  is  the  instantaneous  fluid  velocity,  and  the  integral  is  taken  from  the 
bed  to  the  free  surface.  Theories  presented  below  use  various  approaches  to 
evaluate  the  vertical  distributions  of  the  diffusivity  u,  in  Equation  19  and  u  in 
Equation  22,  and  to  decompose  u  into  mean  and  fluctuating  components  due 
to  currents,  waves,  and  turbulence.  Implicit  in  all  the  following  formulas, 
unless  otherwise  noted,  are  the  following  assumptions: 

a.  The  horizontal  velocity  of  the  sediment  particles  is  equal  to  the 
horizontal  fluid  velocity. 

b.  Flow  accelerations  play  no  explicit  role  in  the  sediment  transport 
processes. 

c.  The  concentration  of  suspended  sediment  is  sufficiently  low  that 
stratification  effects  can  be  neglected. 

As  can  be  seen,  assumptions  for  the  diffusion-based  theories  are  far  fewer 
than  for  the  energetics-based  models,  and  means  to  rationally  extend  the 
theories  are  typically  available,  although  often  analytically  intractable. 


Dally  and  Dean  (1984) 

The  two-layer  model  proposed  by  Dally  and  Dean  for  cross-shore 
suspended  sediment  transport  assumes  that  sediment  motion  occurs  by  two 
mechanisms.  In  the  near-bottom  layer,  each  grain  is  transported  a  small 
distance  before  redepositing  on  the  bed.  Derivation  of  this  "first-order 
contribution,"  in  their  terminology,  assumes  linear  wave  theory  to  predict  the 
mean  horizontal  distance  traveled  by  a  grain  in  one  wave  period.  The  second 
transport  mechanism  is  a  mean  flow  (derived  from  radiation  stress  and 
conservation  of  mass  requirements)  that  operates  in  both  the  near-bed  and 
upper  layers.  The  sum  of  these  sediment  velocity  contributions  is  multiplied 
by  the  sediment  concentration  profile  to  determine  the  vertical  profile  of 
suspended  sediment  transport  using  a  finite-difference  scheme.  The 
concentration  profile  is  modified  from  Rouse  (1937)  to  include  a  contribution 
to  the  shear  velocity  due  to  wave-breaking-induced  turbulence.  Several 
parameters  in  the  theory  are  not  specified  by  Dally  and  Dean,  and  thus 
quantitative  comparison  of  their  theory  with  field  measurements  is  not 
possible. 


Deigaard,  Fredsoe,  and  Hedegaard  (1986) 

Using  the  analogy  of  a  turbulent  bore  to  describe  a  spilling  breaker, 
Deigaard,  Fredsoe,  and  Hedegaard  (1986)  solve  numerically  a  one-equation 
model  describing  the  transport  and  diffusion  of  turbulent  kinetic  energy  in  the 
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surf  zone.  The  total  eddy  viscosity  comprises  contributions  due  to  turbulence 
and  to  the  production  of  energy  in  the  combined  wave-current  motion.  The 
latter  contribution  is  calculated  according  to  the  theory  of  Fredsoe,  Anderson, 
and  Silberg  (1985).  Finally,  the  suspended  sediment  transport  is  calculated 
according  to  Engelund  and  Fredsoe  (1976). 


Semi-Empirical  Models  for  Suspended  Sediment 
Transport 

Given  the  difficulty  in  predicting  suspended  sediment  transport  in  the 
nearshore  from  first  principles,  several  engineering-oriented  approaches  have 
made  extensive  use  of  field  data  from  evolving  beaches  and  from  large-scale 
laboratory  studies  to  parameterize  transport  processes.  The  semi-empirical  or 
schematic  theories  described  below  are  based  on  the  assumption  that  an 
equilibrium  profile  exists  for  any  constant  wave  climate.  Although  that 
assumption  is  currently  the  subject  of  vigorous  study,  the  schematic  approaches 
are  a  practical,  albeit  crude,  means  of  predicting  sediment  transport, 
particularly  when  detailed  fluid-motion  data  are  lacking  and  some  a  priori 
knowledge  of  the  typical  beach  profile  is  available. 


Larson  and  Kraus  (1989) 

Founded  on  considerable  field  and  laboratory  data,  the  SBEACH  model 
developed  by  Larson  and  Kraus  (1989)  (hereafter  referred  to  as  "LK")  extends 
the  energy  dissipation  models  proposed  by  Knebel  and  Dean  (1985)  and 
Kriebel  (1986)  to  predict  total  net  sediment  transport  in  the  surf  zone  and 
foreshore.  These  models  are  designed  to  make  predictions  about  "large-scale 
profile  features  on  intervals  of  tens  of  minutes"  (LK,  p.139).  Thus  they  differ 
fundamentally  in  purpose  from  theories  described  above,  which  predict 
transport  on  a  wave-by-wave  basis.  These  models  make  the  following 
assumptions: 

a.  Beach  profile  change,  and  thus  total  net  sediment  transport,  is  mainly 
governed  by  breaking  of  short-period  waves  having  periods  in  the 
approximate  range  of  3-20  s  (LK,  p.12).  Thus  infragravity  waves  are 
explicitly  assumed  to  have  little  effect  on  sediment  transport  and 
resulting  profile  evolution. 

b.  An  equilibrium  profile  exists  such  that  the  energy  dissipation  per  unit 
volume  is  uniform  across  the  surf  zone  (Dean  1977).  Dean’s  devel¬ 
opment  of  this  hypothesis  assumes  linear  wave  theory,  a  constant  ratio 
of  wave  height  to  water  depth  within  the  surf  zone,  and  linear  and 
parallel  bottom  contours  to  arrive  at  a  characteristic  equilibrium  profile 
shape 
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h  =  Ax 1 

where  A,  the  profile  shape  parameter,  is  a  function  of  particle  size  or  fall 
velocity  and  x  is  distance  offshore.  Departures  from  this  equilibrium  profile 
provide  the  driving  mechanism  for  sediment  transport. 

c.  Waves  are  normally  incident,  and  only  cross-shore  transport  is 
predicted. 

d.  The  direction  of  transport  (i.e.  onshore  or  offshore)  is  determined  by  an 
empirical  criterion  that  differentiates  between  bar  growth  (offshore 
transport),  and  berm  growth  (onshore  transport). 

e.  The  magnitude  of  the  total  net  transport  at  a  point  is  directly 
proportional  to  the  difference  between  the  energy  dissipation  rate  per 
unit  volume  estimated  from  field  measurements  and  the  equilibrium 
dissipation  rate  determined  according  to  Dean’s  (1977)  hypothesis.  The 
transport  formula  proposed  by  LK  (p.  167)  is 


( / )  =  sgn  M  ZZL  - 
WT 


2n  (H0) 


when 

(25) 

and  (/)  =  0  otherwise.  The  sign  or  direction  of  transport  (i.e.,  sgn  )  is 
determined  by  the  quantity  in  curly  brackets,  based  on  an  empirical  criterion 
that  differentiates  between  bar  growth  (offshore  transport  in  the  direction  of 
increasing  coordinate  x)  and  berm  growth  (onshore  transport).  The  quantity 
dh/dx  is  the  bed  slope,  HO  is  the  deepwater  wave  height,  W  is  the  sediment  fall 
velocity,  T  is  the  wave  period  determined  from  the  spectral  frequency  peak  fp, 
M=0.0007  is  an  empirical  constant  and  p,  is  sediment  density.  The  magnitude 
of  transport  is  determined  from  the  quantity  in  the  right-most  square  brackets, 
in  which  E  is  the  energy  dissipation  rate  per  unit  volume  calculated  according 
to  assumptions  discussed  below,  Eeq  is  the  equilibrium  energy  dissipation  rate, 
and  ep  is  a  constant.  The  term  involving  the  beach  slope  eph/dx  is 
incorporated  primarily  to  improve  the  numerical  stability  of  the  model  (LK,  p. 
168).  The  value  of  the  transport  rate  coefficient  K  is  determined  from 
experimental  and  field  measurements,  and  is  discussed  in  more  detail  below. 
The  erosion-only  model  of  Kriebel  and  Dean  (1985)  and  Kriebel  (1986)  is  a 
special  case  of  Equation  24  in  which  transport  is  assumed  to  be  directed 
offshore  only,  and  the  term  involving  the  beach  slope  dh/dx  is  neglected. 
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The  general  form  of  Equation  24  implies  an  interesting  physical  picture  of 
sediment  transport  along  the  cross-shore  profile.  According  to  the  criterion  in 
curly  brackets,  determined  from  quantities  measured  well  outside  the  surf  zone, 
sediment  “is  transported  offshore  or  onshore  along  the  full  length  of  the  active 
profile  at  a  specific  instant  in  time”  (LK,  p.  166).  As  noted  by  LK  (p.  166), 
this  formulation  is  less  likely  to  correctly  predict  the  actual  direction  of 
transport  as  the  beach  approaches  equilibrium,  where  small  departures  from 
equilibrium  may  create  large  qualitative  changes  in  sediment  transport  (i.e., 
change  the  direction  of  transport).  Under  conditions  of  relatively  great  profile 
disequilibrium,  this  deficiency  in  the  model  is  thought  to  be  less  important. 

The  derivation  of  expressions  for  E  and  Eeq  follows  Dean  (1977, 

Appendix  III),  assuming  linear  wave  theory: 

<26> 

h  dxy  8 

where  Hs  is  the  mean  significant  wave  height.  A  useful  form  for  calculating  E 
from  field  measurements  is 


d{H]h) 

Ibc 


(27) 


although  discretization  of  Equation  27  using  less  than  three  relatively  closely 
spaced  measurement  sites  will  yield  predictions  of  transport  at  a  point 
interpolated  between  actual  measurement  points. 


The  expression  for  Eeq  is  (LK,  p.  168) 


where  yb  is  the  breaker  index,  the  ratio  between  wave  height  and  water  depth 
at  breaking.  Equation  28  is  derived  by  substituting  ybh  for  Hs  in  Equation  26 
and  differentiating  with  respect  to  x  to  obtain 


Eeqdx= 


which  can  be  integrated  to  give  Equation  23. 
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The  profile  shape  parameter  A  is  then 


A  = 


2 

24  Eeq 


PS7  T2 


(30) 


which  can  in  turn  be  solved  for  Eeq. 

Dean  (1987)  proposed  an  empirical  dimensional  formula  for  determining  A 
as  a  function  of  settling  velocity  W 

A  =  0.067W044  '  (31) 

in  which  A  has  units  of  meters1/3  when  the  settling  velocity  W  is  expressed  in 
cm  s'1.  The  fit  to  the  data  used  by  Dean  (1987,  p.  11)  appears  good  for 
D50  >  1  mm,  and  less  good  for  smaller  grains  likely  to  be  transported  in 
suspension.  It  is  not  clear  whether  LK  use  Equation  31  to  determine  A  or  the 
earlier  relationship  proposed  by  Moore  (1982). 


SBEACH  parameters 

Various  parameters  for  SBEACH  have  been  determined  for  laboratory  and 
field  application.  The  particular  choice  of  wave  height,  for  example,  depends 
on  the  intended  application.  This  section  summarizes  and  discusses  parameters 
and  measurements  relevant  to  field  application  of  SBEACH. 

The  recommended  value  of  the  transport  rate  coefficient  K  appearing  in 
Equation  24  was  determined  for  both  laboratory  and  field  situations  by 
minimizing  the  sum  of  squares  of  deviations  between  measured  and  predicted 
profile  elevations.  Seven  erosional  (bar-forming)  profiles  generated  in  large 
wave  laboratory  facilities  were  used  to  determine  an  optimal  value  for  K  of 
1.6  10-6  m4  AT1;  in  the  field,  a  value  of  K=Q.l  10'6  m4  V1  minimized  the 
discrepancy  between  measured  and  calculated  profiles.  Moore  (1982) 
determined  a  value  of  K=  2.2  lO-6  m4  V1  using  laboratory  and  field  data  in  a 
similar  model  (no  explicit  term  proportional  to  dhldx).  Intercomparison 
between  SBEACH  and  the  model  of  Kriebel  (1986)  produced  an  optimal 
K=8.1  10'6  m4  AT1  for  the  Kriebel  model.  Because  it  is  not  possible  to  relate 
the  predicted  transport  rate  to  any  physically  measureable  quantity,  finding  an 
optimal  calibration  procedure  for  K  is  difficult.  Finally,  the  breaker  index  yb 
is  set  equal  to  unity  for  field  application  (LK,  p.  207).  Note  that  yb  differs 
from  y,  the  ratio  of  wave  height  to  water  depth  across  the  surf  zone  (Thornton 
and  Guza  1983). 
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Ackers  and  White  (1973) 


The  semi-empirical  sediment  transport  model  of  Ackers  and  White  (1973) 
was  developed  for  steady,  uniform  open-channel  flows,  and  is  considered  here 
because  it  is  widely  used  in  a  variety  of  applications.  The  model  describes 
total  sediment  transport  in  terms  of  three  dimensionless  groups  related  to  the 
size,  mobility,  and  transport  of  grains  having  diameters  greater  than  0.04  mm. 
Nearly  1,000  flume  experiments  using  well-sorted  sediment  in  steady,  uniform 
flows  having  depths  up  to  0.4  m  form  the  empirical  basis  for  the  model.  A 
single  test  of  the  model  under  a  quasi-unidirectional  flow  regime  was  per¬ 
formed  during  one  tidal  cycle  in  an  estuary  (p.  2053): 

...  theoretical  predictions  appear  to  exaggerate  the  net 
transport  in  a  landward  direction.  This  is  probably  a 
result  of  the  asymmetrical  nature  of  the  ebb/flood  cycle 
and  there  is  clearly  a  need  for  further  investigation  when 
the  flow  is  time-dependent. 

Use  of  the  Ackers  and  White  formulation  is  thus  not  recommended  for 
nearshore  oscillatory  flows.  The  following  sections  outline  the  procedure  for 
computing  sediment  transport  of  fine  grains  in  suspension.  A  dimensionless 
grain  diameter  D„  mobility  parameter  F„  and  transport  parameter  G,  are 
calculated  as  follows: 

The  dimensionless  grain  diameter  D,  is 


D  ■  D 


8{  PrP) 


ptr 


(32) 


where  D  is  the  medium  grain  size  and  v  is  the  kinematic  viscosity.  Ancillary 
parameters  n,  Ag,  Cg,  and  m  are  empirically  derived  functions  of  D .: 


n  =  1  -  0.56  log 


(33) 


_  i 

A  =  0.23  D, 7  +  0.14 

g  * 


(34) 


logC?  =  2.86 log D,  -  ( logD,  )2  -  3.53 


(35) 
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m  =  9.66  D'1  +  1.34 


(36) 


Taking  n~  1  for  fine  suspended  sediments,  the  mobility  parameter  F,  is 
simply  the  square  root  of  the  Shield’s  parameter  0 


F  = 


[S0(pf-p]’ 


=  eT 


(37) 


and  the  transport  parameter  G.  is  then 


G 

* 


(38) 


Ackers  and  White  formulate  the  sediment  flux  in  terms  of  the  sediment 
concentration  in  parts  per  million  by  weight  of  fluid  flux,  and  thus  the 
immersed  sediment  transport  rate  is 


6  P  J  P  r 

(  V 

F 

_ 1  -i 

P  ,Du 

P,  1 

<A*  , 

[pk  J 

(p  uh ) 


(39) 


where  u  is  the  depth-averaged  velocity,  and  the  various  factors  of  h,  p  and  p5 
have  been  left  unsimplified  for  clarity.  For  field  application,  the  shear  velocity 
u,  in  Equation  39  might  be  computed  using  Equation  4  and  a  suitable  choice 
for  cf . 
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Chapter  2  Theory 


Table  1 

Transport  Formulae 


Bowen  (1980)1 


Bailard  (198 1)2 


(r)=pc/-^ 


Roelvink  and  Stive  (1989)  (cross-shore  transport  only) 


where  AT,  =p  [  exp  (h/H^ )  -  1  ]"'  f 


Larson  and  Kraus  (1989)  (cross-shore  transport  only) 


,  .  1  {H0)  3  2n(H0)\  ep  dh 

(p'“p,s  [ 


when  E  >  E  -  £■_  —  ,  and  zero  otherwise. 
K  dx 


Ackers  and  White  (1973) 


(i) =  10"6 


p  -p  F 
—L-LC  —  -  1 


P  sDu 


(P  uh) 


Bagnold  (1963)  is  functionally  equivalent,  but  not  elaborated  in  the  same  detail. 
Humiston  (1993)  follows  Bailard  (1981). 
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Appendix  A 
Notation 


The  notation  adopted  herein  attempts  to  use  commonly  accepted  symbols 
whenever  possible.  Coordinate  systems  are  all  right-handed,  but  in  some  cases 
the  ^-coordinate  is  parallel  to  the  sloping  planar  bed  while  in  others  the  x-y 
coordinate  plane  is  horizontal.  Except  where  noted,  angle  brackets  are  used  to 
denote  a  temporal  or  ensemble  average;  overbars  indicate  a  spatial  average; 
tildes  indicate  an  oscillatory  component;  and  primes  indicate  a  fluctuating 
component.  The  subscripts  b  and  0  refer  to  quantities  evaluated  at  the  bed, 
although  the  location  of  the  nominal  bed  elevation  differs  in  various  theories. 
The  symbol  [0]  denotes  a  dimensionless  quantity. 


c 

Wave  phase  velocity 

[LT1] 

cf 

Friction  factor  or  bed  drag  coefficient 

[0] 

c* 

Group  velocity 

[LT1] 

Cr 

Coefficient  in  Equation  21 

[0] 

dhldx 

Bed  slope 

[0] 

eb 

Bed-load  efficiency  factor 

[0] 

es 

Suspended  load  efficiency  factor 

[0] 

f 

Wave  frequency 

[T1] 

fP 

Peak  wave  frequency 

rr1] 

L 

Wave  friction  factor 

[0] 

g 

Gravitational  acceleration 

[L  T-2] 

h 

Water  depth 

[L] 

i 

Immersed  weight  sediment  transport  rate  per  unit  width 

[ML1  T1] 
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*0 

Immersed  weight  sediment  transport  rate  in  8  direction 

[ML'1  T1] 

Unit  vector  in  ^-direction  (offshore) 

[0] 

£ 

Wave  number  (  =  2  n  L'1) 

[L-1] 

*, 

Turbulent  kinetic  energy  per  unit  mass 

[L2  T2] 

m 

Exponent  defined  in  Equation  36 

[0] 

n 

Exponent  defined  in  Equation  33 

[0] 

s 

Ratio  of  grain  and  fluid  densities 

[0] 

t 

Time 

[T] 

u 

Water  velocity  in  x-direction 

[LT1] 

u. 

Friction  or  shear  velocity 

[LT1] 

Um 

Maximum  orbital  wave  velocity 

[LT1] 

«e 

Current  velocity  in  0-direction 

[LT1] 

u 

Instantaneous  fluid  velocity 

[LT1] 

V 

Water  velocity  in  y-direction 

[LT1] 

w 

Water  velocity  in  z-direction 

[LT1] 

X 

Horizontal  coordinate  (+offshore) 

[L] 

y 

Horizontal  coordinate 

[L] 

z 

Vertical  coordinate  (+up) 

[L] 

A 

Profile  shape  parameter 

[L1/3] 

Constant  defined  in  Equation  34 

[0] 

B 

Coefficient  in  Equation  13 

[0] 

Constant  defined  in  Equation  35 

[0] 

D 

Mean  sediment  grain  diameter 

[L] 

D. 

Dimensionless  sediment  grain  diameter 

[0] 

E 

Energy  dissipation  rate  per  unit  volume 

[M  L'1  T-3] 

A2 
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Equilibrium  energy  dissipation  rate  per  unit  volume 

[M  L-1  T3] 

F, 

Dimensionless  mobility  parameter 

[0] 

G. 

Dimensionless  transport  parameter 

[0] 

H 

Wave  height 

[L] 

Ho 

Deepwater  wave  height 

[L] 

Hs 

Significant  wave  height 

[L] 

Hrns 

Root-mean-square  wave  height 

[L] 

L 

Wave  length 

[L] 

K 

Coefficient  in  Equation  24 

[L3  M1  T2] 

Kb 

Coefficient  in  Equation  3 

[L1  T2] 

K, 

Coefficient  in  Equation  6 

[L3  T-2] 

k2 

Coefficient  in  Equation  6 

[L4  T  ] 

M 

Coefficient  in  Equation  24 

[0], 

N 

Volume  concentration 

[0] 

Nbed 

Volume  concentration  of  randomly  packed  particles 

[0] 

S 

Energy  slope  in  unidirectional  flow 

[0] 

R 

Rouse  number 

[0] 

T 

Wave  period 

[T] 

ux 

Perturbation  velocity 

[LT1] 

w 

Sediment  fall  velocity 

[LT1] 

a 

Wave  angle  (deviation  from  normal) 

[radians] 

P 

Bed  slope 

[0] 

P* 

Coefficient  in  Equation  1 1 

[0] 

Energy  dissipation  rate 

[M  T-3] 

ep 

Coefficient  in  Equation  24 

[M  T3] 
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T 

Ratio  of  wave  height  to  water  depth 

[0] 

yb 

Ratio  of  wave  height  to  water  depth  at  breaking 

[0] 

K 

von  Karman’s  constant  (=  0.4) 

[0] 

0 

Shield’s  parameter 

[0] 

P 

Fluid  density 

[M  L-3] 

P, 

Sediment  density 

[M  L-3] 

o 

Radial  wave  frequency 

[T1] 

% 

Bed  shear  stress 

[M  L'1  T'2] 

V 

Kinematic  viscosity 

[L2  T1] 

Sediment  diffusivity 

[L2  T1] 

CO 

Wave  or  stream  power 

[M  T'3] 

A4 
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